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Abstract 


This  note  is  a  description  of  the  subprograms  DOP , 
which  calculates  the  Doppler  velocity  of  an  earth-bound 
observer,  JULDA  which  calculates  the  Julian  day-number, 
MOVE  which  calculates  precession,  COORD  which  performs 
coordinate  transformations,  and  the  GRM  series  of  sub¬ 
programs  which  deal  with  Doppler  velocities  in  terms  of 
a  standard  galactic  rotation  model.  These  subprograms 
are  written  in  basic  FORTRAN  and  should  be  useable  on  a 
wide  variety  of  computers. 


Accepted  for  the  Air  Force 

Franklin  C.  Hudson 

Chief,  Lincoln  Laboratory  Office 


SUBROUTINE  OOP  (RAHRs.RAmin.RASEC  .DUEC  *DMIN.DSF'~  .NYR.fif  ..Y, 

2  NHUT  »  NMUT  *  NSUT  *  ALAT  »  OLONG  *  ELEV »  XLST * VSUN . VMOri . VOBS » V 1 > 

C 

C  OOP  CALCULATES  THE  VELOCITY  COMPONENT  OF  THE  OBSERVER  *1 1 TH  RESPECT 
C  TO  THE  LOCAL  STANDARO  OF  REST  AS  PROJECTED  ONTO  A  LINF  SPECIFIED  RY  THr  *IGHT 
C  ASCENSION  AND  DECLINATION  (RAHKS.  RAMIN.  R  ASEC  .  JDCG.  DM  IN.  L'StC)  EPOCH  OF 
C  DATE.  FOR  A  TIME  SPECIFIED  AS  FOLLOWS*  NYR  =  LAST  TWO  DIGITS  OF  THE  YFaR 
C  (FOR  19aX  A.D.).  NDAY  =  DAY  NUMBER  ( G  -4T )  »  NHUT*  NMUT*  NSUT  =  HRS «  MJn,  SEC 
C  (GMT).  THE  LOCATION  OF  THE  OBSERVER  IS  SPECIFIED  BY  THE  LATITUDE  (ALAT)* 

C  L0N3ITUJE  (OLONG)  (GEODETIC)  (IN  DEGREES)  AND  ELEVATION  (FLEV)  ( T»>  METERS) 

C  ABOVE  MEAN  SEA  LEVEL.  THE  SUBROUTINE  OUTPUTS  THE  LOCAL  mEan  SlDr°EAL  T IMF 
C  (XLST  I,*  DAYS)*  THE  COMPONENT  uF  THE  SUN'S  MOTIOI  aITH  RESPECT  TO  THE  LOCAL 
C  ST  A  JDARu  OF  REST  AS  PROJECTED  ONTO  THE  LINE  OF  SIGHT  TO  THE  SOURCE  (VSUN  IN 
C  KM/SEC)  as  well  as  THE  TOTAL  VELOCITY  COMPONENT  VI  (KM/SFC).  POSITIVF 
C  VELOCITY  CORRESPONDS  TO  INCREASING  DISTANCE  BETWEEN  SOURCE  AMD  OBSERVFP. 

c 

C  THIS  VERSION  OF  OOP  TAKtS  INTO  ACCOUNT  COMPONENTS  OF  THF  0  SFMVFR'S 
C  MOTION  JUE  TO  THt  ROTATION  OF  THE  EARTH.  THE  REVOLUTION  OF  THF  EARTH-MCOD 
C  6ARYCENTER  A30UT  THt  SUN.  AND  THE  MOTION  OF  THE  EARTH'S  CENTFf!  ABOUT  THE 
C  EARTH-MOON  BARYCENTER*  THE  PERTURBATIONS  OF  THE  EARTH'S  OPMT  DUF  10  THE 
C  PLANETS  ARE  NEGLECTED.  THE  ABSOLUTE  PRECISION  OF  THIS  VERSION  OF  OOP  IS 
C  ABOUT  0.004  KM/SEC.  bUT  SINCE  THE  DOMINANT  ERROR  TERM  IS  SLOWLY  VARYING.  THF 
C  RELATIVE  ERROR  WILL  8E  CONSIDERABLY  LESS  FOR  TIMES  UP  TO  A  -TEK  OP  SO. 

C 

C  REFERENCES*  MCRAE*  D.  A.*  wESTERHOUT .  G..  TABLE  FOR  THE  REDUCTION  OF 
C  VELOCITIES  TO  THE  LOCAL  STANDARD  OF  REST.  THF  OBSFRVAT  ;RY . 

C  LUND.  SWEDEN.  1966. 

C  SMART.  W.  M..  TEXT-BOOK  ON  SPHERICAL  ASTRONOMY.  CAmorieGE 

C  UNIV.  PRESS*  1962. 

C  THE  AMERICAN  EPHEMERIS  AND  NAUTICAL  ALMANAC 

C  THE  oUPPLEMENT  TO  THE  ABOVE 

C 

C  VERSION  OF  JUNt  1969 

C 
C 
C 
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c 

C  THE  FOLLOWING  CALCULATIONS  DEAL  .-(ITh  THr  SUN’S  *’  >TION  XITh  ASPECT  TO  THF- 
C  uOCAL  STANDARD  OF  REST. 

L  THE  VlLuCITY  OF  Tml  SUN  W ITH  REoPcCT  TO  THE  LOCAL  STANDARD  OF  RFST  IS  THE 
C  CONVENT IOivAL  VALUE  OF  20.0  KV/SEC  TOrtAR;  RA  =  IB  HRS »  DEC  =  30  DEG  (  1900). 
AAA=18. 0*3. 1413926535/12.0 
DD=3o. 0*3. 14 15926535/ 180.0 
C  <10VE  PROCESSES  THIS  DIRECTION  TO  DATE 

CALL  MOVE  (190U>1900+NYRrl»NJAY*AAA»DD*  DEL  A  * OELDO > DC )  MOVF 

AAA=AAA+DELA 

DD=UD+DELDD 

i.  THIS  VELOCITY  IS  CONVERTED  TO  CARTESIAN  COMPONENTS 
X0=20. 0*COSF ( A AA ) *COSF ( DD ) 

YO=20.0*SIiNF(AAA)*CGSF(DO) 

ZO=2j.O*SlNF (Du> 

C  RA1  IS  THE  RIGHT  ASCENSION  (REVS=)AYS> 

R  Al  =  ( RAHRS+R  AM  IN/60  •  0+RASEC/3600 . 0  )  /  24  ,  O' 

C  RA  IS  Tiit  RIGHT  ASCENSION  (RADIANS) 
i<A=2.0*3. 1415926535*RA1 
C  DEC  IS  THFI  DECLINATION  (RADIANS) 

JEC  =  3. 141 5926535* (DUtG  +  SIGNF(U ^IN/60.0  +  DSEC/360 0.0 i DDrcO  ) / 180.0 
C  CC.  CS»  AMD  S  ARE  THE  DIRECTION  COSINES  CORRESPONDING  TO  RA  A\G  D^C 
CC=CuSF ( DEC ) *COSF ( RA ) 

CS=CoSF (DEC ) *SINF (RA ) 

S=SINF (DEC) 

C  VSUJ  IS  Tml  PROJECTION  ONTO  THE  LINE  OF  SIGHT  TO  TmE  STAR  OF  THE  SUN’S 
C  ,«OTIOM  WITH  RtSPECT  TO  THE  LOCAL  STANDARD  OF  REST  (KM/SEC) 
VSUN=-XO*CC-YO*CS-ZO*S 

c 

C  COORDI'iATlS  OF  THE  ObSERVtR»  LATITUDE  (RADIANS) i  AND  LONGIY  JnE  (RFVS=D/.YS) 
LAT=ALAT*3.14l392o535/180.0 
•MLON>j=OLONG/360 . 0 

c 
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o  r. 


the  fullowing  calculat ions  deal  ;th  the  ti 'E 

C  THE  EPOCH  IS  1900  JANUARY  O.b  UT  =  JULIAN  DAY  241S'i20.0 

C  JU  lb  T lit  TIME  FROM  THE  tPOCH  TO  JAN  0.0  OF  THE  CUKRFNT  YEAR  (JAYS) 

UU= ( JULDA ( 1900+NYK ) -241S020 ) -0 .b  J  iLP A 

C  TU  lb  DJ  CONVERTED  TO  JULIAN  CENTURIES 
TU=DJ/Jbb2b.O 

C  UTJA  IS  THE  GMT  FROM  JAN  0.0  TO  THE  PRESENT  ( DAY  j ) 

UTDA=NDAY-HvlHUT/24.0  +  NMUT/1440.0  +  NSUT/864  00.0 

C  SMD  (SMALL  D)  IS  THL  TIME  FROM  T HE  EPOCH  TO  THE  PR  SENT  (DAYS) 

SMD=JU+U  TOA 

C  T  IS  SMJ  CONVERTED  To  JULIAN  CENTURIES 
T— SM^/3b32S.O 

C  START  IS  THE  GREENWICH  MEAN  SIDEREAL  TIME  ON  JAN  0.0  (JAYS) 

C  (THE  EXTRA  129.1794  SECS  CORRESPONDS  TO  THE  0.7  CENTURY  SuHTR ACTFD  FROM  TU. 

C  THE  PRECISION  IS  THEREBY  IMPROVED.) 

ST«K  T=  (b.0  +  ,3».0/b0.0+ (4b.  83o+129. 1794  + A64  0  184.5m  2*  (  T'J-0 . 7 ) +0 . 0^29 
2  *TU**2)/3bOU.O)/24.0 

C  Cl  IS  THE  CONVERSION  FACTOR  FROM  SOLAR  TIME  TO  SIDEREAL  TIME 
Cl=(23.0+Sb.0/b0.U+4.090S4/3b00.0>/24.0 

C  GST  IS  THE  GREENWICH  MEAN  SIDEREAL  TIME  ( PAYS > 
oST-bTAkT+UTDA/tl 

C  XLST  IS  THl.  LOCAL  MEaN  SIDEREAL  TIME  (F(<OM  JAN  0)  (DAYS) 

XLST=GST-WLONG 

XLST=XLST-IFIX(XLST) 

C 
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C  T HL  FOLoOjv  1NG  CALCULATIONS  JFAL  /ITH  THE  OBSERVER'S  'OTIOI.  /.TTH 

C  RESPECT  TO  THE  tARTH'S  CENTER. 

C  REDUCTION  Of-  GEODETIC  LATITUDE  TO  GEOCENTRIC  LAT1T  *DE  ( ARCSFC  ONDS ) 
JLAT=-<11.0*b0.04D2.743U>*SlNF<2.0*CAT)  +  1.1633*'>I^F(4.(l*Cf'T) 

2  -0.002o*SInM6.0*CAT) 

C  CONVERT  CAT  TO  GEOCENTRIC  LATITUDE  (RADIANS) 
CAT=CAT4ULAT*3.l41b92oS35/3b00.U/180.0 
C  Rho  IS  THE  RADIUS  VECTOR  FROM  THE  EARTH'S  CENTER  To  THE  CuSPRVEP  C  •'■ETE  .  S> 
R  -o37blb0.0*  (0.998o27073  +  U.n016  76438*C0SF(2.0*CAT)-n.0l'00(J3SlO 
2  *CO_>F(4.0*CAT)+O.OOOUOOOU4«COSF<6.0*CAT) )4FLIV 
C  AND  VRHO  IS  THE  CORRESPONDING  oIRCULAR  VELOCITY  ( Mt.TFRS/S  IHEr.  EAL  HAY) 
VRn0=2.0*3.14lS92bb3S*Rri0 
C  CONVERTED  TO  K ILOMETLRS/SEC 

VRHO=VRHO/24.0E3/3bOU.O*Cl 

C  VOES  IS  THE  PROJECTION  ONTO  THc.  LINE  OF  SIGHT  TO  THE  STAP  CP  THE  VELOCITY 
C  OF  THE  OBSERVER  WITH  RESPECT  TO  THE  EARTH'S  CENTER  (KM/SEC) 

VObS=VRHU*COSF ( C AT ) *COSF ( DEC ) *SINF ( 2 . 0*3 . 1 4 1 59?b53S* ( XLST-R A 1 ) ) 

C 
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C  THE  FOLLOWING  CALCULATIONS  '">£ AL  wITH  THE  EARTH’S  ORBIT  ftRijtiT  THE  SUN 
C  AM  IS  THE  MEAN  ANOMALY  (OF  THE  EARTH'S  ORBIT)  (RADIANS) 

AS  =  <  05c).  R75»3*U.  9BSt>UU,2b70*SMn-U.  000  lS0*T**2-0. '10000  3*t**:m 
2  *3.1415926535/180.0 

C  PI  IS  TnE  MEAN  LO.xGITUCt  OF  PERIHELION  (RADIANS) 

PIr< 101. 220B3*J.0000470684*SMn*0.000453*T**2+0. 30000 3*T**3) 

2  *3.1415926535/180.0 

C  E  IS  THil  ECCENTRICITY  OF  THE  ORBIT  (DIMENSIONLESS) 

E=u. U 1675 104-0. 00004 180* T-0. 000000 126* T* *2 
C  A I  IS  TnE  MEAN  OBLIQUITY  OF  THE  ECLIPTIC  (RADIANS) 

Air (23.452294-u.0130125*T-0. 00000164*T**240.000H00503*T**3) 

2  *3.1415926535/180.0 

C  VS  IS  THE  TRUE  ANOMALY  (APPROXIMATE  FORMULA)  (RADIANS) 

C  (EQUATION  OF  THE  CENTER) 

VS=A>i-M2. 0*E-0.25*E**3)  *SINF  (  AM)+i.25*E**2*SINF(2.0*A^)* 

2  13.U/12.0*E**3*SINF(3.0*AM)  ♦. . . 

C  XLAM  IS  THE  TRUE  LONGITUDE  OF  THE  EARTH  AS  SEEN  FROM  THE  SUN  (RADIANS) 
XLAMrp I ♦VS 

C  ALAM  IS  THE  TRUE  LONGITUDE  OF  THE  SUN  AS  SEEN  FROM  THE  EARTH  (RADIANS) 
ALAM=XLAM43. 1415926535 

C  BETA  IS  THE  LATITUDE  OF  THE  STAR  (RADIANS) 

C  ALONG  IS  THE  LONGITUDE  OF  THE  STAR  (RADIANS) 

CALL  COORD  ( 0 . 0 » 0 . 0 » -3. 1415926535/2 . 0 » 3. 1415926535/2 . 0-A I . 

2  R A  »  DEC  »  ALONG  t BETA )  COORD 

c  aa  is  the  semi-major  axis  of  the  earth’s  orbit  <km) 

AArl49598500.0 

C  an  is  the  mean  angular  rate  of  the  earth  about  THE  SUN  (RATI A-iS/DAY) 

ANr2. 0*3. 1415926535/365.2564 

C  HOP  IS  H/P  FROM  SMART  =  THE  COMPONENT  OF  THE  EARTH’S  VELOCITY  PERPENDICULAR 
C  TO  THE  RADIUS  VECTOR  (Km/DAY) 

HOP=AN*AA/SQRTF  ( 1 .0-.,* *2) 

C  CONVERTED  TO  KM/SEC 
HCP=HOP/86400.0 

C  V  IS  THE  PROJECTION  ONTO  THE  LINE  OF  SIGHT  TO  THE  STAR  OF  THr  VELOCITY 
C  OF  THE  lARTH-MOON  BARYCEMTER  WITH  RESPECT  TO  THE  SUN  (KM/c^fC) 

Vr-n^P*C0SF (BETA)* (SINF ( ALAM-ALONG)-E*SINF (PI-ALOUG) ) 

C 


5 


c 

C  THE  FOLLOWING  CALCULATIONS  DEAL  .  ITH  THE  BOON'S  OR;  IT  AROU’  D  THE 
C  tARTh-MOON  BARYCENTER 

C  OMGA  (OrlEGA)  IS  THE  LONGITUDE  UF  THE  MEAN  ASCENDING  NODE  OF  the  LUNAR  ORBIT 
C  (DEGREES) 

OMGA  =259.  1832 75-0.  (j529539222*SWD+0. 002078*T**2+0 .  0  00" 02*T**3 
C  OMGAK  Is  OMGA  IN  RADIANS 

OMGAR=OMGA*3.1415926535/180.0 

C  AMON  Is  OMGA  PLUS  THE  MEAN  LUNAR  LONGITUDE  OF  THE  MOON  (DEgRFFS) 

C  (SHOULD  BE  13 . 176396526b ) 

AMON=270. 434164+13. 176396527  *SMD-0 . 001 1 33*T**?+0 . 00000 1R*T**3 
C  GAMP  (GAMMA-PRIME)  IS  OMGA  PLUS  THE  LUNAR  LONGITUDt  OF  LUNAR  PERIGEE  (DEGREFS 
GAMP=334. 329556+0. 1114040803*SMD-0. 01 0325*T**2-n. 00 00 12*T**3 
C  PIM  IS  THE  MEAN  LUNAR  LONGITUDE  OF  LUNAR  PERIGEE  (TO  RADIANS) 

P I M=(GAMP-OMGA) *3.1415926535/180.0 
C  EM  IS  THE  ECCENTRICITY  OF  THE  LUNAR  ORBIT 
EM=0. 054900489 

C  OLAMM  Is  THE  MEAN  LUNAR  LONGITUDE  OF  THE  MOON  (TO  RADIANS) 

OLAMM= ( AMON-OMGA  >*3.1415926535/180 . 0 
C  AIM  IS  THE  INCLINATION  OF  THE  LUNAR  ORBIT  TO  THE  ECLIPTIC  (RADIANS) 

A I M=b. 1453964*3.1415926535/180.0 
C  AMM  IS  THE  APPROXIMATE  MEAN  ANOMALY  (RADIANS) 

C  (IT  IS  APPROXIMATE  bECAUSE  PIM  SHOULD  bE  THE  TRUE  RATHER  THAU  THE  MEAN  LUNAR 
C  LONGITUJE  OF  LUNAR  PERIGEE) 

AMf»i=OLAMM-PIM 

C  VSM  IS  THE  TRUE  ANOMALY  ( APPROXIMATE  FORMULA)  (RADIANS) 

C  (EQUATION  OF  THE  CENTER) 

VS  -,=AMv+(2.0*EM-0.25*EM**3)  *SINF  ( AMM )  +1  »25*EM**2*SINF  ( 2 . 0* AMM ) 

2  +13.0/12.0*EM**3*SINF(3.0*AMM)  +... 

C  ALAMM  IS  THE  TRUE  LUNAR  LONGITUDE  OF  THE  MOON  (PADIANS) 

ALAMm=PIM+VSM 

C  ANM  IS  THE  MEAN  ANGULAR  RATE  OF  THE  LUNAR  ROTATION  ( R  AO  I  A'.S/O  A  Y ) 

ANM=2. 0*3. 14 15926535/27. 32 lbol 

C  A AM  IS  THE  SEMI-MAJOR  AXIS  OF  THE  LUNAR  OBRIT  (KILOMETERS) 

AAm=oO. 2665*6378 .388 

C  BETAM  IS  THE  LUNAR  LATITUDE  OF  THE  STAR  (RADIANS) 

C  ALG M  IS  THE  LUNAR  LONGITUDE  OF  THE  STAR  (RADIANS) 

CALL  COORD  ( OMGAR » 0 . 0  * OMGAR-3. 1415926535/2 . 0 » 3 . 1415926535/2 . 0- AIM. 

2  ALONG  »bET A  *  ALGM. BETAM)  COORD 

C  HOPM  IS  H/P  FROM  SMART  =  THE  COMPONENT  OF  THE  LUNAR  VELOCITY  PERPENDICULAR 
C  TO  ThE  RADIUS  VECTOR  (KM/DAY) 

HOPM=ANM*AAM/SGRTF ( 1 • 0-EM**2 ) 

C  CONVERTED  TO  KM/SEC 
HOPM=HOPM/86400.0 

C  VMON  IS  THE  PROJECTION  ONTO  THE  LINE  OF  SIGHT  TO  THE  STAR  OF  THE  VELOCITY 
C  OF  THE  EARTH’ S  CENTER  WITH  RESPECT  TO  THE  EARTH-MOON  BARYCENTER  (KM/SEC) 

C  (THE  81.30  IS  THE  RATIO  OF  THE  EARTH’S  MASS  TO  THE  MOON’S  ha<;s) 

VM0N=-H0PM/8l .30*C0SF (BETAM) * (SINF ( ALAMM-ALGM)-+M*SINF (PlM-ALOM) ) 

C 
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Vl=V+VSUN+VMQN+VOl  S 
RETURN 


THlb  PROGRAM  OMITS  THE  PLANETARY  PERTURBATIONS  ON  THE  EAPTF  K  QPniT.  THr 
AMOUNT  TO  ABOUT  O.OUJ.  KM/SEC  AND  ARE  THOUGHT  TO  BE  THE  LARGEST  CONTRI BUT  I 
TO  THE  ERROR  IN  THE  VELOCITY. 


Li'ji  ) 
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it  c 


function  julda (nyr) 

C  THIS  FUNCTION  COMPUTES  THE  JULIAN  DAY  NUM8EP  AT  12  HRS  UT  ON  JANUARY  0 

C  OF  THE  YEAR  NYR  (GREGORIAN  CALENDAR).  JULDA  IS  AN  INTEGER  BFCAUSE  OF  THIS 
C  DEFINITION.  FOR  EXAMPLE*  JULDA  =  2439856  FOR  NYR  =  1968. 

C 

UYRM1=NYR-1 
IC-NYRM1/ 100 

JULDA=1721425+365*NYRMl+NYRMl/4-IC+IC/4 

RETURN 

END 
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o  r>  o  r  o  r  r>  o  r;  r*  r>  o  c  r  r>  r>  r>  o  o  r>  r"  r>  n  r  r>  r  0000000000 


SUBROUTINE  COORD  ( AO » BO »  AP»  BP.  '  1  *  81  >  A2.B  2) 

THIS  oUoROUTlNt  CONVERTS  THE  LONGITUDE-LIKE  (Al)  AND  LA  T II  U'JE-L  IKF  (Rl) 
COORDINATES  OF  A  POINT  ON  A  SPHEKF.  INTO  THE  CORRESPONDING  COORDINATES  (A?. 

62)  IN  A  DIFFERENT  COORDINATE  SYSTEM  THAT  IS  SPECIFIED  BY  INF  COORDINATES 
OF  ITS  ORIGIN  (AO.  BO)  AND  ITS  NORTH  POLE  (AP.  BP)  IN  THE  ORIGINAL  COOl DINATf 
SYSTEM.  THE  RANGE  OF  A2  WILL  BE  cROM  -PI  TO  PI. 

ALL  AKG JMEN  rs  ARE  IN  RADIANS. 

EXAMPLE.,  OF  USE 

PI  =  3. 1415926S3:> 

PI02  =  PI/2.0 

EXAMPLE  I — TO  CALCULATE  AZIMUTH  AND  ELEVATION  FR  )M  HOUR  ANGLE  AND  UECLINAT ION 
CALL  COORD  (PI »P 102-LATITUDE. 0 .0 .LATITUDE* HOUR  ANGLE . DE CLI NATI on, 

2  AZI  lUTH.ELEVATION) 

f HZ  4  IF  AZIMUTH  IS  DESIRED  IN  THE  RANGE  0  TO  PI  SET 
AZIMUTm  =  AZIMUTH  ♦  (PI  -  SIGNFCPI .AZIMUTH)) 

lXAMPLc.  II  —  TO  calculate  hour  angle  AND  DECLINATION  from  AZIMUTH  and 
ELEVATIoN 

CALL  COORD  (PI rPI02-LATITUDE» 0,0. LATITUDE. AZIMUTH.ELEVATIC 
2  HOUR  ANGLl.DLCLINATION) 

EXAMPLE  III — TO  CALCULATE  LI  AND  Dl  FROM  RIGHT  ASCLNSION  AND  DECLINATION 
(EPOCH  1900.0) 

AP  =  (12. 0+40. 0/60. 0)*PI/12.0 
(I.E.  12  HOURS  40  MINUTES  CONVERTED  TO  RADIANS) 

BP  =  26 «  0*P I /ldO  •  0 

AO  =  ( 18.0+40.0/o0.0) *PI/12.0 

BO  =  0.0 

(REFER  TO  KRAUS.  P.»  RADIO  ASTRONOMY.  MCGRAW  HILL.  NEW  YORK.  1966.  BUT  FOR 
FURTHER  REFINEMENTS.  SEE  ALSO  ALLEN.  C.  W. >  ASTROPnYSIC AL  QUANTITIES. 

ATHLONE  PRESS.  LONDON.  1963.) 

CALL  COORD  ( AO.BO. AP.BP.RIGHT  ASCENSION.DECLI JATION.LI >t I ) 
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000000000000000000000000000000000 


EXA  *1PLc.  IV — To  CALCULATE  RIGHT  ASCENSION  AND  DECLINATION  (EPOCH  1900.0)  FROM 
LI  AND  Bl 

IN  GENERAL*  whenever  we  KNOW  THE  FORWARD  TRANSFORMATION  (EXAMPLE  III  ABOVE) 

WE  MAY  DO  THE  REVERSE  TRANSFORMATION  WITH  AT  MOST  TWO  EXTRA  PRELIMINARY  CALLS 
TO  COOKu  TO  CALCULATE  THE  COORDINATES  IN  SYSTEM  2  OF  THE  POLF  AND  ORIGIN  IN 
SYSTEM  1.  BUT  OFTEN  IT  IS  POSSIBLE  TO  GET  THESE  NEEDED  COORDINATES  BY 
INSPECTION.  FOR  EXAMPLE*  BP  WILL  REMAIN  THE  SAME  FOR  THE  FORWARD  AND  REVERSE 

transformations,  for  this  example  we  see  by  inspection  that 

APP  =  6. 0*P 1/12.0 
BPP  =  28 • 0*PI / 180 • 0 

(THE  SECOND  P  REPRESENTS  PRIME.)  AND  WE  MAY  CALCULATE  AOP  AriD  BOP  FROM 
CALL  COORD  < AO  *  BO  * AP , BP  *  0 . 0  *  0 . 0  * AOP *  BOP ) 

WHERE  THE  AO*  ETC.  ARE  FROM  EXAMPLE  III.  THEN  THE  ACTUAL  CONVERSION  IS 
CALL  COORD  ( AOP *  BOP  * APP * BPP *  LI  * BI *  RIGHT  ASCENSION*DECLINATION) 

tXAMPLE  V — TO  CALCULATE  LI  I  AND  BII  FROM  RIGHT  ASCENSION  AND  DECLINATION 
(EPOCH  1950.0) 

AP  =  (12.0+49.0/oU.U)*PI/12.0 
BP  =  27.4*PI/180.0 
AO  =  (17. 0+42. 4/60. 0)*PI/12.0 
BO  =  -(28. 0+55. 0/60. 0)*PI/l80.0 

CAlL  COORD  ( AO*BO* AP,BP*RIGHT  ASCENSION* DECLI NATION  *  LI  I . BI I ) 

EXAMPLl  VI — TO  CALCULATE  RIGHT  ASCENSION  AND  DECLINATION  (EPOCH  1950.0)  F70M 
LI  I  AND  BII 

first  calculate  app  and  bpp  from 

CALL  COORD  ( AO  *  BO  * AP , BP  *  0 . 0  * PI02 * APP * OPP > 

THEN  CALCULATE  AOP  AND  BOP  FROM 

CALL  COORD  ( AO  *  BO  * AP *  BP  *  0 . 0  *  0 . 0  * AOP *  BOP) 

WHERE  THE  AO*  ETC.  ARE  FROM  EXAMPLE  V.  THEN  THE  ACTUAL  CONVFRSION  IS 
CALL  COORD  ( AOP *  BOP  * APP * BPP * L 1 1  * BI I  *  RIGHT  ASCENSION.DECLIMATION) 
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tXAMPLE  VII  —  TO  CALCULATE  (ECLIPTIC)  LATITUDE  AND  LONGITU  r  p*OM  RIGHT 
ASCEri^IvN  AND  DECLINATION 

lPS  i  ->  THE  OBLIGUITY  OF  THE  ECLIPTIC  WHICH  IS  ABOUT  23. '*43  DEGREES*  BUT  IT 
C  OEPE  .oS  ON  THE  EPOCH.  SEE  THE  AMERICAN  EPHE^ERIS  Ain  .  AUT 1 0  AL  ALMANAC. 

C  t-PS=^3.443*PI/IbO.C 

C  CALL  COOKD  ( 0 . 0 . 0 . 0  * -P 102 . P 1 02-tPS . R I GHT  ASCFNSI OU*  DECLI  NATION. 

C  2  LATITUDE. LONGITUDE) 

C 

C 

C  THE  '  uTAT ION  USES  S  OR  C  FOR  SINt  OR  COSINE  OF  THE  CDRRESPCN  ING  VAR  I  At  LF  > 

C  FOR  EXAMPLE.  SBO  =  SINF(HO).  ETC. 

C 

C  NOTE  that  THE  INPUT  PARAMETERS  ARE  PARTIALLY  REDUNDANT.  FOR  EXAMPLE »  IF 
C  APr  BP«  AND  AO  ARl  SPECIFIED.  THEN  THERE  ARE  ONLY  TWO  DISCRETE  VALUES 
C  POSSIBLE  FOR  BO  (EXCEPT  FOR  A  FEW  DEGENERATE  SPECIAL  CASFS).  SFE  BELOl  FOR 
C  <vHAT  TO  DO  IF  IT  IS  NECESSARY  TO  PRECALCULATE  AO  ADD  BO. 

C 

C  IF.  INSTEAD  OF  AO  AND  bO.  THE  LONGITUDE  OF  THE  ASCENDING  NODF  IS  KNOWN  I' 

C  BOTH  THt  OLD  ( AN1  )  AND  NEW  ( AN2 5  COORDINATE  SYSTEMS.  THEN  AO  AND  HO  '■‘A'i  ME 
C  CALCULATED  bY  A  PRELIMINARY  CALL  TO  COORD 
C  CALL  COORD  ( 0 . 0 . 0 . 0 » AN1-AP . BP » -AN2 . U . 0 . AO . BO ) 

C  THEN  THIS  AO  AND  BO  MAY  BE  USED  FOR  A  SERIES  OF  ORDINARY  CALLS  TO  COOPD  AS 
C  DESCRIbED  AbOVE. 

C 

C  IF  AP.  bP.  AND  AO  ARE  KNOWN.  THEN  THE  TWO  POSSIBLE  VALUES  OF  RG  may  BE 
C  CALCULATED  FROM 

C  SbO=(SBP  $  2.0*CbP**2*CAPAO*SGRTF ( I .0+CAPAO**2) ) / 

C  2  (SdP**2+(CBP*CAPA0)**2) 

C  WHERE  CAPAO  =  CCSF(AP-AO)  AND  THE  OTHER  NOTATION  IS  EXPL'-I'F  AFOVfc  .  Ai  D 
C  WHERE  THE  5  IS  TO  BE  REPLACED  bY  +  AND  -. 

C 

C  IF  AP.  bP*  AND  bO  ARL  KNOWN.  THEN  THE  TWO  POSSIBLE  VALJFS  OF  AO  v A Y  BF 
C  CALCULATED  FROM 

C  CAPAu=(1.0-SbO*SdP)/<C80*C8P) 

C  BOTH  ANGLES  WITH  THIS  COSINE  ARE  POSSIBLE. 

C 
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SBO=SINF (BO) 

CBO=COSF(BO) 

SbP=SINF(BP) 

CHP=COSF(bP) 

SB1=SINF(B1 ) 

Col=COSF(Bl) 

C 

SB2=bBP*SBl+CbP*CBl*C0SF (AP-A1  ) 

D2=AblNF(SB2) 

C  (NOTE  bO  IS  NOT  NEEUtU  TO  CALCULATE  R2) 

CB2=C0SF (B2) 

C 

SAA=SINF ( AP-A1 ) *Cbl/CB2 
CAA=(SHl-SB2+SbP)/(CB2+CBP) 

C 

CdB=bBO/CBP 
bBB=slNF (AP-AO) *CbO 
C 

TA202=( 1 .0-CAA*CQ8-SAA+SBB)/ (SAA*CBB-CAA*SBB) 
A2=2.0*ATANF(TA202) 

C 

KETUPN 

END 
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SUBROUTINE  MOVE  < NYR I *NYRF *  MO »  .DA *RA *C *  UELR  *  HELD  *  f)C  ) 

C 

C  MOVE  CALCULATES  THl  CORRECT  I OH  (DELHI  IN  RIGHT  ASCENSION  (  .A)  AND  THE 
C  CORRECTION  <UELU)  IN  DECLINATION  (D)  (ALL  IN  RADIANS)  TO  ME  ADJFD  TO  ThE 
C  mean  COORDINATES  FOR  EPOCH  NYRI  (  .G.  19bn)  TO  GIVE  THE  APPARENT  POSITIONS 
C  OF  A  DATE  SPECIFIED  BY  THf  YEAR  ( NYRF *  E.3.  1968)*  MONTH  CO*  1  TO  12),  AMD 
C  DAY  (NOA).  IF  THE  DAY-N JMBER  IS  KNOWN*  USE  IT  FOR  NQA  AND  SET  -  1. 

C  MOVt  Al.jO  CALCULATES  THE  t'OUATIOM  OF  THE  FOUINOXLS  (DC.  IN  *  I  MUTES  OF  TI'F) 

C  WHICH  AaY  BE  ADUED  TO  THE  MEAN  SIDEREAL  TIME  TO  6IVC  THE  APPARENT  SIDEi  EAL 
C  TIME  (AtNA-169).  UELR  AND  DELD  CONTAIN  CORRECTIONS  FOR  PRECESSION*  ANNUAL 
C  ABERRATION*  AND  SOME  TERMS  OF  NUTATION.  IF  RA  AND  U  ARE  FOR  THE  MEAN  tPOCH 
U  (I.E.  HALFWAY  BETWEEN  NYRI  AND  NYRF)  THEN  THE  PRECISION  OF  NFLR  AND  OELU  IS 
C  ABOUT  2  ARCSECONDS  (SEE  NEGLECTED  TERMS  IN  ESE-44).  IF  RA  AND  D  ARf  EITHER 
C  OF  THL  -NO  POINTS  OF  THE  INTERVAL*  THEN  THE  PRECISION  MAY  BE  SOMEWHAT  WORSE. 
C  AENA  =  THE  AMERICAN  EPHEMER IS  AND  NAUTICAL  ALMANAC  (THE  BLUE  BOOK). 

C  ESE  =  THE  EXPLANATORY  SUPPLEMENT  TO  ABOVE  (THE  GREEN  BOOK). 

L 

SNO= jINF (D ) 

CSD=COSF(D) 

TND=SNO/CSO 

c 

uSR=oOSF (RA) 

SNR=SINF(RA) 

C 

C  AL  IS  AN  APPROXIMATE  DAY  NUMBER  (I.E.  THF  NUMBER  OF  DAYS  SINCE  JANUARY  0 
C  OF  THE  YEAR  NYRF). 

AL=3j*(MO-1)+NDm 

C 

C  TO  IS  THE  TIME  FROM  1900  TO  NYRI  (CENTURIES) 

TO=Fl_OATF(,n|YRI-1900)  / 100.0 

C  T  IS  THE  TIME  FROM  NYRI  TO  DATE  <  lYRF  *  MO*  NDA)  (CENTURIES) 

C  (365. 2421988  IS  THE  NUKbER  OF  EPHFMERIS  CAYS  IN  A  TROPICAL  YEAR) 

T=(Fi_OATF  (NYRF -NYRI )  +AL/365. 2421988 )/100.0 
C  2ETA0  IS  A  PRECESSIONAL  ANGLE  FRO-'  ESE-29  (ARCSECONDS) 

ZET Au=< 2304. 250+ 1. 396*T0 ) *T+0.302*T **2+0.0 lfi*T* *3 
C  DITTO  FuR  Z 

Z=ZETAO+0.791*T**2 
C  AND  THETA 

THET A=( 2004. 682-0. 853*TO ) *T-0. 426* T**2-0. 042* T**3 
C  AM  AND  AN  ARE  THE  M  AND  N  PRECESSIONAL  NUMBERS  (SEE  AEUA-bu*  474)  (RAOlAiS) 
AM=UtTAO  +  Z)  *4. 64813681  IE-6 
AN  =  T  )ETA*4 , 848 13681 IE-6 
C 
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ALAM  IS  AN  APPROXIMATE  MEAN  LONGITUDE  FOR  THE  SU *  ( AENA-SO )  (RADIANS) 
A|_AM=:(  0.985647*  AL+278. 5)  *0.0174532925 
SNL=SINF( ALAM) 

CSL=COSF( ALAM) 

C  DELR  IS  THE  ANNUAL  ABERRATION  TERM  In  RA  (RADIANS)  (ESE-47.4P) 

C  (0.9174S051  =  COSIOBLIGUITT  OF  ECLIPTIC)) 

C  ( -9 . 924 1360 5E-5  =  K  =  20.47  ARCSECONDS  =  CONSTANT  OF  ABERRATION  (ESE-4P)) 
DELR=-9.92413605E-5*(SNL*SNR+0.91745051*CSL*CSR)/CSD 
C  PLUS  PRECESSION  TERMS  (SEE  AENA-50  AND  ESE-38) 

2  +AM+AN*SNR*TND 

C  DELD  IS  DITTO  ABOVE  IN  DECLINATION 

DELD=-9. 9241 3605E-5* (SNL*CSR*SND-0. 91 745051 *CSL*SNR*SND 
C  (0.39784993  =  SIN ( OBLIQUITY  OF  ECLIPTIC)) 

2  +0.39784993*CSL*CSD)  +AN*CSR 
C 

C  THE  FOLLOWING  CALCULATES  THE  NUTATION  (APPROXIMATELY)  (ESL-41,45) 

C  OMEGA  IS  THE  ANGLE  OF  THE  FIRST  TERM  OF  NUTATION  (ESE-44)  (APPROXIMATE 
C  FORMULA)  (DEGREES) 

0MEGA=259. 183275-1934 . 142* ( TO+T ) 

C  ARG  IS  OMEGA  CONVERTED  TO  RADIANS 
ARG=OMEGA*0. 01 74532925 

C  DLONG  IS  THE  NUTATION  IN  LONGITUDE  (DELTA-PSI)  (RADIANS) 
0L0NG=-8.3597E-5*SINF( ARG) 

C  QOBLQ  IS  THE  NUTATION  IN  OBLIQUITY  (DELTA-EPSILON)  (RADIANS ) 

DOBLQ=  4.4678E-5*C0SF( ARG) 

C 

C  ADD  NUTATION  IN  RA  INTO  DELR  (ESE-43) 

UELR=DELR+DLONG* (0.91745051  +0.39784993  *SNR*TND) -CSP*TND*DOPLO 
C  AND  DEC. 

DELD=DELD+0. 39784993  *CSR*DLONG+SNR*OORLQ 
C  DC  IS  THE  EQUATION  OF  THE  EUUINOXES  (MINUTES  OF  TIME)  (ESE-43) 

DC  =  Di_ONG*21  0 . 264  lb9 

RETURN 

END 
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SUBROUTINE  GRWIl  (ELII.VR.  RS.^C.PSI) 

DIMENSION  RS(2) »PSI (2) 

C  j JJKOUT iNE  GRM1  (GALACTIC  ROTATION  MODEL  -  1)  ACCEPTS  THE  GALACTIC  LONGITUDE 
C  (ELI  I  Im  OtGREES)  AND  THE  DOPPLER  VELOCITY  WITH  RESPECT  To  THc  LOCAL  STANDARD 
C  OF  REST  (VR  IN  KM/SEC)  ANL  RETURNS  THE  DISTANCE  FROM  THE  SuN  (RS) *  THE 
C  JISTANCc  FROM  THE  GALACTIC  CENTER  (RC)  (BOTH  IN  KPC).  AND  THE  GALACTOCFNTRIC 
C  LONGITUDE  (PSI  IN  DEGREES).  THE  RANGE  OF  PSI  WILL  QE  -IPO  TO  +  1R0  DEGREES. 

C  FOR  MANY  VALUES  OF  ELII  AND  VR .  THERE  WILL  BE  A  DISTANCE  AMBIGUITY  AND  THE 
C  1WO  POSSIBLE  DISTANCES  WILL  BE  RETURNED  IN  RS ( 1 )  AND  RS(2)  ANC  THr  ASSOCIATED 
C  ANGLES  In  PSI(l)  AND  PSI (2) •  IF  RS ( 1 )  IS  NEGATIVF  OR  ZERO.  THEN  THERE  IS 
C  ONLY  ONc.  bOLUT  I  ON  AND  IT  IS  IN  RS(2)  AND  PSI  ( 2 )  •  IF  ROTH  PS(1)  AND  RS(2) 

C  ARE  NlohTIVE  OK  ZERO.  THEN  THERE  ARE  NO  SOLUTIONS  TO  THE  EQUATIONS.  I.E.  THF 
C  GIVEN  VR  IS  IMPOSSIBLE  FOR  THIS  ELII. 

C  SEE  SU3.\0UTINE  GR -13  FOR  FURTHER  DETAILS. 

C 

PI =3. 1415926535 
TPI=2.0*PI 

C  ELI IR  IS  ELII  CONCERTED  TO  RADIANS 
ELIIr(=ELII*PI/180.0 
C  SLII  IS  S IN ( ELI  I ) 

SL1I=SINF(ELIIR> 

C  RO  IS  T ,,E  SUN'S  DISTANCE  FROM  THE  GALACTIC  CENTER  (KPC) 

RO=lU.O 

C  OMGAO  lb  THE  ANGULAR  ROTATION  VELOCITY  AT  THE  SUN  (KM/SEC/K°C) 

OMGAu=25.0 

C  OMGA  IS  THE  ANGULAR  ROTATION  VELOCITY  AT  THE  SOURCE  (KM/SEC/kPC) 

OMGA=VR/ (RO*SLII ) +OMGAO 
C  GRM3  CALCULATES  RC  (KPC) 

CALL  GRM3  (OMGA.  RC)  GRM3 

C  H  IS  THE  LENGTH  OF  THE  PERPENDICULAR  DROPPED  FROM  THE  GALACTIC  CENTER  TO  THE 
C  LINc.  UF  SIGHT  (KPC) 
h=RO*SLII 

C  GO  TO  SO  IF  THERE  IS  NO  SOLUTION 
IF  (ABSF(H)-RC)  20.20.SO 

C  JELT  IS  THE  ANGLE  BETWEEN  THE  SUN  AND  THE  GALACTIC  CENTER  AS  SEEN  FROM  T  <r 
C  SOURCE  (RADIANS) 

20  DELT=ASINF(H/RC) 

C  PSI1  AND  PSI2  ARE  THE  TWO  SOLUTIONS  FOR  PSI  (RADIANS) 

PSI 1-ELI 1R-UELT 
PSI2=ELIIR+DELT+PI 

C  CONVERT  THESE  TO  DEGREES.  SET  THE  RANGE  TO  -180  TO  180.  AND  PUT  IN  PSI 
PSI ( 1)=PSI1* 180.0/PI -360.0W IF  I X( PSI 1/TPI+SIGNF( 0.5 »PSI1) ) 

PSI (<i)=PSI2*lb0.0/PI-360.0*IFIX(PSI2/TPI+SIGNIF (0.5. PSI 2) ) 

C  THEN  SET  RS ( 1 )  AND  RS(2)  (KPC) 

RS ( 1 )=-RC *SI NF (PSI 1) /SLII 
RS(2)=-RC*SINF(PSI2)/SLII 
RETURN 
C 

C  GO  TO  SO  FOR  NO  SOLUTIONS  TO  THE  EQUATIONS  (A(3SF(H/RC)  .GT.  1) 

50  HS(1)=0.0 
RS ( 2 ) =0 . 0 
PSI ( 1 ) =0 • 0 
PSI (2)=0.0 
RETURN 
END 
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sUdR JUT  INC  GRM2  (ELI1#RS»  VR#RC»PSI) 

C  SUBROUTINE  GR'ia  (GALACTIC  ROTATION  MODEL  -  2)  ACCEPTS  THF  <  AI.ACTIC  LONr-IT'lDF 
C  C  El_  1 1  I  -I  DEGREES)  AND  ThE  DISTANCE  FROM  THE  SUN  <RS  IN  KPC )  AND  RFTURNS  THE 
C  DOPPLER  VELOCITY  WITH  RESPECT  TO  THE  LOCAL  STANDARD  OF  REST  (VR  IN  KM/SEC)* 

C  THE  DISTANCE  FRON  THE  GALACTIC  CENTER  <RC  IN  KPC)«  AND  THE  GALACTOCENTE IC 
C  LONuI  TUJc.  (PSI  IN  DEGREES). 

C  SEE  SUBROUTINE  GKM4  FOR  FURTHER  DETAILS. 

C 

PI=i. 141592653S 

C  ELI IR  la  ELI I  CONVERTED  TO  RADIANS 
ELI IH=ELI I *P 1/160.0 

C  RO  IS  THE  SUN'S  DISTANCE  F HOP  THE  GALACTIC  CENTER  (KPC) 

RU=10.0 

C  THEN  RC  IS  JUSI 

RC=SJRTF  <RS**2+R0**2-2.U*RS*R0*C0SF(ELIIR) ) 

C  JMGA  (PROP  GRK4)  IS  THE  ANGULAR  ROTATION  VELOCITY  ;«T  THE  SC  'JPCfc  ( KM/SEC/K^C ) 
CALL  GRN4  <RC»  OMGA )  GKN4 

C  uMGAO  la  THE  ANGULAR  ROTATION  VELOCITY  AT  THE  SUN  <KH/SEC/kpC) 

UMG A  d— 25 . 0 

C  AND  THEN  VR  IS  JUST 

VR=(OMGA-OMGAO) *RO*SINF(ELIIR) 

C  a  IS  THE  ANGLE  BETWEEN  THE  SUN  AND  THE  SOURCE  AS  SEEN  FROM  ThF.  GALACTIC 
C  CENTER  (RADIANS) 

C  SNA  IS  aIN(A)  (DIMENSIONLESS) 

SNA=RS*SINF(LLIIK)/RC 
C  CSA  IS  COS ( A )  (DIMENSIONLESS) 

CSA=(R0**2+RC**2-RS**2) /<2.U*R0*RC> 

C  THEN  A  IS  JUST 

A=2.0*ATANF(SNA/(1.0+CSA) ) 

C  AND  PSI  (CONVERTED  Ju  DEGREES)  IS  JUST 
PSI=-A*loO.G/F I 
RETURN 
END 
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SUBROUTINE  GKM3  (UMGA.  RC ) 

jURKOU  r  1 NE  URM3  (bMLACTIC  ROTATION  MODEL  -  3)  ACCEPTS  THE  ASHLAR  ROTATION 
VELOCITY  (OMGA  IN  KM/StC/KPC)  AND  RETURNS  THE  DISTANCE  FROM  THE  GALACTIC 
CENTER  (KC  IN  KPC ) .  BASED  ON  THE  MODEL  OF  THE  GALAXY  SUGGESTED  BY  SCHMIDT 
( 19ob ) •  b£E  SUBROUTINE  GRM4  FOR  FURTHER  DETAILS. 

IF  GMbA  ,LE.  0  THEN  RC  =  0 
IF  (oMGA)  lu.lO»2U 
10  RC=0.0 

return 

UNFORTU JATeLY  THERE  iS  NO  WAY  TO  INVERT  SCHMIDT'S  EQUATIONS  TO  OBTAIN  j<C 
DIRECTLY  FROM  OMGA.  THE  FOLLOWS  4  ITERATION  PROCEDURE  SEEMS  TO  WORK  AS  WELL 
AS  A.gY .  IT  IS  uASED  ON  THE  APPROXIMATION  THAT  THE  CIRCULAR  VFLOCITY 
(=OMGA*RC>  IS  A  SLOWLY  VARYING  FUNCTION  OF  RC.  THl  NUMBEk  BOO  IS  A  FIoST 
GUESS  AT  THIS  VELOCITY. 

20  RC  =  <2.iO.O/OMGA 

Du  3u  1  =  1.5  I 

CALL  GRM4  (RC.  OMGAT)  GRM4 

(NOTE  TiimT  H  GOES  TO  1  FOR  RC  MUCH  LESS  THAN  10.  AND  TO  2/3  FOR  LARGE  RC ) 

B=  (<2.0+10,0/(10*  O+RC  )  )/3.0 
RC=RC+RC*B* (OMGAT-OMGA) /OMGA 

30  CONTINUE  I 

RETUnN 
END 
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SUijR^UT  INE  GRM4  (RC*  OMGA ) 

C  SUBROUT INE  GRM4  (GALACTIC  ROTATION  MODEL  -  4)  ACCEPTS  THE  (  T STANCE  FROM  THE 
C  GALACTIC  CENTER  (RC  IN  KRC)  AND  RETURNS  THE  ANGULAR  ROTATION  VELOCITY  (O^GA 
C  IN  KM/SEC/KPC)  .  BASED  ON  THE  MODEL  OF  THE  GALAXY  SUGGESTED  BY  MAARTEN 
C  SCHMIDT  *  CHAP,  IN  GALACTIC  STRUCTURE*  ED.  BY  A.  BLAAUW  A  in  ,M,  SCHMIDT* 

C  UNIV.  OF  CHICAGO  PRESS*  1965. 

C 

C  IF  RC  .LE.  0  THEN  OMGA  =  0 
IF  (,<C>  10  *  10  *  2u 
10  OMGA.1U.U 
return 
c 

20  IF  (kC-10.0)  40*30*50 
30  OMGA-25 . U 
RE  TURN 


the  following  calculations  are  for  rc  .lt.  10  kpc 

VC  IS  The  CIRCULAR  VELOCITY  (KM/SEC) 

40  VC=SjRTF (30000. 0/RC+10120.2*RC-41.722*RC**3) 
THEN  OMGA  IS  JUST 
OMGA=VC/RC 
return 

the  following  calculations  are  for  RC  .GT.  10  KPC 

50  VC=S jRTF( 85161 1 .6/RC-214858S. 1/RC**2 ) 

omga=vc/rc 

return 

END 
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13  ABSTRACT 


This  note  is  a  description  of  the  subprograms  DOP,  which  calculates  the 
Doppler  velocity  of  an  earth-bound  observer,  JULDA  which  calculates  the 
Julian  day-number,  MOVE  which  calculates  precession,  COORD  which  per¬ 
forms  coordinate  transformations,  and  the  GRM  series  of  subprograms  which 
deal  with  Doppler  velocities  in  terms  of  a  standard  galactic  rotation  model. 
These  subprograms  are  written  in  basic  FORTRAN  and  should  be  useable  on 
a  wide  variety  of  computers. 
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